A search for incipient lattice instabilities in MgB2 by anelastic spectroscopy 



F.Cordero/ R. Cantelli,^ G. Giunchi,^ and S. Ceresara^ 
^ CNR, Istituto di Acustica "O.M. Corbino", Via del Fosso del Cavaliere 100, 
1-00133 Roma, Raly and Unita INFM-Romal, P.le A. Mom 2,1-00185 Roma, Raly 
^ Universitd di Roma "La Sapienza", Dipartimento di Fisica, and 
Unita INFM-Roma 1, P.le A. Mora 2, 1-00185 Roma, Raly 
3 EDISON S.p.A, Foro Buonaparte 31, 1-20121 Milano, Raly 

We measured the complex dynamic elastic modulus of MgB2 between 1.3 and 650 K at frequencies 
included between 5 and 70 kHz. The main purpose of this work is to search for possible lattice 
instabilities, which could contribute to enhance the electron-phonon coupling. The Young's modulus 
indeed presents an anomalous softening on cooling below 400 K, and hysteresis between cooling 
and heating. Such anomalies, however, can be accounted for by intense relaxation processes with 
maxima between 50 and 150 K, whose nature is not yet clarified, and by geometric defects generated 
by the anisotropic thermal expansion between the grains during thermal cycling. Additional intense 
relaxation processes are observed above room temperature, and their possible origin is discussed. 
The infiuence of the highly anharmonic in-plane vibration modes of the B atoms on the elastic and 
anelastic properties of MgB2 is discussed in detail. 
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The discovery of superconductivity with Tc = 39 K 
in MgB2 is arousing great interest, since this intermetal- 
Uc compound seems to be a simple BCS superconductor, 
without any of the other mechanisms for high-Tc super- 
conductivity that are hypothesized for the p.uprate su- 
perconductors. Electronic band calculationsElo indicate 
a predominant hole conduction in the two-dimensional u 
band from the spxPy orbitals of B, with an extremely 
large deformation potential for the B bond stretching 
mode, which dominates the electron-phonon coupling. 
The most recent first priBciple calculations of the elec- 
tronic bands and phononsa indicate that the zone-center 
optical mode with symmetry, involving vibrations 
of the B atoms within the ah plane, has very flat anhar- 
monic potential, and can account for the observed Tc with 
a calculated electron-phonon coupling constant A ~ 1. 
An experimental indication of the anharmonicity of that 
mode comes from the high value of its Griineisen parame- 
ter, evaluated from the shift of the corresponding Raman 
band with pressureJj Still, it is important to ascertain 
whether there are other anharmonic or unstable modes 
that are not easily seen by phonon DOS measurementsci 
and other spectroscopies, but can contribute to the large 
electron-phonon coupling. In many cases of BCS su- 
perconductivity with high Tc, the strongly anharmonic 
phonons causing the enhanced electron-phonon coupling 
also produce a structural instability, which in general pro- 
foundly affects the acoustic properties of the material. 
The best known examples are the A15 superconductors,Q 
where a clear correlation exists between the occurrence of 
a transformation from cubic to tetragonal, often assimi- 
lated to a martensitic transformation, and a high value 
of Tc. It should be emphasized that the occurrence and 
the entity of the structural instability in the A15 com- 
pounds is rather sensitive to the material preparation, 
and the easiest method for revealing it was the obser- 
vation of softening of the elastic moduli on cooling be- 
low room temperature.Q Also the superconducjtijig per- 
ovskite BKBO presents structural instabilitiesjj'tl while 
ultrasonic absorption peaks in BPBO are interpreted in 
terms of strong phonon anharmonicityQ 

In the case of MgB2 a possible unstable mode is the 
Big mode of the out-of plane vibrations of the B atoms, 
even though the first jprinciple calculations indicate a stiff 
and harmonic mode.H In fact, it has been observed that 
MgB2 is close to the condition of buckled B planes in a 
phase diagram of the AB2 compounds versus charge den- 
sity and B-B distance.cl In addition, when the hole density 
in the B planes is progressively reduced by partially sub- 
stituting Mg with Al, the c lattice parameter collapses 
in correspondence with the disappearance of supercon- 
ductivity, suggesting that MgB2 is close to a. structural 
instability at slightly higher electron densityJj 

The anelastic spectroscopy is very sensitive to struc- 
tural instabilities, and therefore we measured the com- 
plex dynamic Young's modulus of MgB2 between 1.3 and 
650 K, searching for possible signs of lattiee instabili- 
ties not yet detected by neutron diffractioncHl3 and other 



spectroscopies. 

The sample was prepared by reacting Mg and B in a 
sealed stainless steel tube, lined with Nb foil, at 950 °C 
for 2 houps, Further preparative details are reported 
elsewhere. Ilil A bar 39 x 6 x 1 mm'^ was cut by spark ero- 
sion. The mass density was p — 2.27 g/cvp? , 86% of the 
theoretical one. The superconducting transition, mea- 
sured by resistivity, was centered at 39.5 K, with a width 
of 1 K. The X-ray diffraction spectra showed the presence 
of a residual minority phase of metallic Mg, presumably 
at the grain boundaries. The sample was suspended with 
thin thermocouple wires and electrostatically excited on 
its fiexural modes, whose resonant frequencies are given 
by /„ = anh/f^E'/p, where E{lj,T)^ E' + iE" is 
the complex dynamic Young's modulus,E3 I and h the 
sample length and thickness, and a„ numerical constants 
depending on the vibration mode. The 1st, 3rd and 5th 
modes could be excited, at uj/2Tr = 5.7, 38 and 73 kHz 
respectively at 1 K. The value of the Young's modulus at 
room temperature, not corrected for the sample poros- 
ity and the Mg traces, is E = 167 GPa. The imaginary 
part is related to the elastic energy loss coefficient Q^^ — 
E" /E' , which was measured from the width of the reso- 
nance curves. 

Figure 1 shows the anelastic spectrum between 1 and 
620 K, measured by exciting the first fiexural mode. A 
first result is the absence of anomalous softening be- 
low 150 K, so excluding a scenario like that of the 
A15 compounds, where the strong phonon anharmonic- 
ity is due to a mode which becomes unstable slightly 
above Tc. Still, anomalous softening and hysteresis be- 
tween cooling and heating are observed between 400 and 
150 K. Such anomalies, however, should not be straight- 
forwardly attributed to a structural instability, since 
microplastic phenomena connected with the anisotropic 
thermal expansion of the hexagonal crystallites are ex- 
pected. Neutron diffractior£3 experiments have shown a 
large anisotropy of the thermal expansion and compress- 
ibility, both being about twice larger along the c axis. 
These plastic phenomena have already been observed, 
in other hexagonal polycrystals, like Zn, Co and Cd,E£l 
whose thermal expansivity is sufficiently anisotropic and 
the critical stress ctc for plastic deformation sufficiently 
low to yield the formation of dislocations during ther- 
mal cycling. The increase of the density of geometrical 
defects results in a reduction of the modulus, and the 
competition between the defect generation and annihila- 
tion, the latter especially at high temperature, explains 
the hysteresis between cooling and heating above 130 K. 
Below that temperature the thermal expansion becomes 
so small that the stress generated at the grain boundaries 
is below (Jc and no further modification in the geometrical 
defects occurs. We do not know what type of geometri- 
cal defects accommodates the strains generated during 
thermal cycling in MgB2, but they are not exclusively 
accommodated at the Mg impurities between the grains; 
in fact, after heating the sample up to 690 K in vacuum, 
about 8% of Mg was lost, the absorption component at 
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80 K was reduced, but the hysteresis of the modulus be- 
low room temperature was unchanged. Also the elastic 
energy loss is partially due to the motion of the geomet- 
rical defects, at least above 150 K, since it has a compo- 
nent which increases with the cooling/heating rate, and 
exhibits hysteresis, as shown by the open symbols in Fig. 
1. The plastic phenomena can account for the irregu- 
larities and hysteresis of E (T) between 150 and 400 K 
but probably not completely for the plateau of E' (T) 
between these temperatures. This plateau could also 
have another origin; indeed, a small softening of some 
modes appears in the phonon DOS measured by neutron 
scatteringEl around 325 K. We shall see however, that 
also the relaxation processes observed below room tem- 
perature in the elastic energy loss might cause a similar 
effect. 

The intense absorption peak at 475-|K is due to a 
thermally activated relaxation process,t3 namely it can 
be described as a contribution 6E = —A/ (1 + iajr) to 
E {u!,T), where the relaxation time follows the Arrhenius 
law T = To exp (W/T) . The absorption peak and mod- 
ulus dispersion of this type of processes shift to higher 
temperature with increasing frequency. Due to the high 
damping, the higher frequency modes could be followed 
only up to 440 K, but exhibited such a shift, and could 
be interpolated by the expression 



SE = -A/[l + {iujTf] 



(1) 



where the parameter a — 0.65 reproduces the broaden- 
ing of the spectrum of the energy barriers (a = 1 in the 
monodispersive case), W = 0.89 eV and tq = 10~^^ s. 
The dotted lines in Fig. 1 are the corresponding contri- 
butions to the absorption and to E' , plus the tail of an- 
other process at higher temperature, probably connected 
with grain-boundary motion or with the loss of Mg. The 
value of To is typical of point defect motion, as could be 
the migration of Mg atoms and vacancies, but, due to 
the bad quality of the higher frequency data, we cannot 
exclude that tq is actually longer and therefore connected 
with the motion of geometrical defects. 

In view of the occurrence of plastic phenomena dur- 
ing thermal cycling above ~ 130 K, it is likely that the 
absorption contains an important contribution from the 
relaxation of extended lattice defects. Some contribution 
may also be expected from Mg at the grain boundaries, 
since deformedJVIg presents an extremely broad absorp- 
tion maximurrO reminiscent of the curves found here 
below room temperature; nevertheless, the contribution 
of metallic Mg is weighted with its molar fraction, i.e. 
few percents, and ccmsidering that the intensity of the 
peak in deformed MgO never exceeds 2 x 10^^, it cannot 
account for elastic energy loss curves of Fig. 2. In addi- 
tion, at Tc it is evident a change of the slope of both the 
imaginary and real parts of E (T) , which must be con- 
nected with processes occurring in the superconducting 
bulk of MgB2. Another unusual feature of the anelastic 
spectrum below 200 K is that the curve simultaneously 



measured at higher frequency, instead of being simply 
shifted to higher temperature, as for any thermally acti- 
vated relaxation process, is about a factor 1.32 more in- 
tense than that at lower frequency. Finally, it should be 
noted that the spectrum below 100 K is perfectly repro- 
ducible during various cooling and heating cycles; con- 
sidering the extreme sensitivity of the dislocation peaks 
to the sample state, and the occurrence of plastic phe- 
nomena during thermal cycling, it is unlikely that the 
low temperature spectrum is exclusively connected with 
geometrical defects. 

It is tempting to relate the low temperature anelas- 
tic spectrum with the highly anharmonic opti- 
cal mode which is thought responsible for the strong 
electron-phonon coupling. Indeed, a relaxation mech- 
anism involving the soft optical modes responsible for 
the structural phase transfpanations in perovskites has 
been proposed by Barrett, E3 and later adopted in or- 
der to explain the ultxasonic absorption in the per- 
ovskite superconductorQ BaRbj _a;Bia;03 and even in 
some cuprate superconductors.E3 The mjechanism is anal- 
ogous to the Akhiezer-type absorption,t3 but instead of 
the acoustic phonons, it is based on the existence of 
optical phonons with a large Griineisen constant, 7 = 
— d(lna;o) /de, namely whose frequency loq strongly de- 
pends on a strain e. The flexural vibration of the sam- 
ple causes a non homogeneous uniaxial strain e along the 
sample length and modulates the energies of the phonons 
according to their Griineisen parameters. The phonons 
system goes out of thermal equilibrium if different modes 
have different values of 7, and the equilibrium is es- 
tablished with an effective phonon relaxation time Tcff. 
The change of the population is reflected in the dynamic 
elastic modulus, which acquires an out of phase imag- 
inary component and a temperature dependent soften- 
ing. It is appropriate to discuss this mechanism, since 
the Griineisen constant of the F,2g optical mode has beep, 
found to have a very large value at room temperature, u 
~ i^Tna — being the lattice constant, while 

normal values range between 1 and 2. Barrett solved the 
case in which the system is schematized with an optical 
mode with frequency and large 7op, and the acoustic 
modes with normal 7ac; since we are considering uniax- 
ial strain s = da /a instead of volume change, we define 
7 = 3 (7op — 7ac) — 7, where 7ac ~ 1.5 is assumed. The 
contribution of the piionon relaxation to the dynamic 
Young's modulus E isEJ 
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where vq is the cell volume, Cac and Cop are the specific 
heat of the acoustic and optic modes per molecule and 
Tac and Top are their relaxation times. The elastic energy 
loss coefficient is, for luTcs ^ 1, 



3 



E' 



Q-=^^ ^TCopUJT^s = A (T) wTeff ; (4) 



where 



a 



op 



voE 
has 



the 



sual 



expression 



fee [{hujQ/2T) / 8mh{hujo/2T)] withB hujo/kB = 830 K, 
E = 167 GPa and vq = 29 x 10^24 ^^^3. result- 
ing relaxation strength A is 6.6 x 10^^ at 100 K and 
6.5 X 10""^ at 300 K, but has to be multiplied by the 
factor LOTcff, which is very small. In fact, the phonon re- 
laxation rates and therefore T~f^ should exceed 10^^ s~^; 
Tac can be deduced from the lattice thermal conductivity 
K through the collision formula k — -ICacW^Tac- Assum- 
ing K ~ 0.1 W/(cm K) foundliZl above 100 K, an average 
so und velocity somewhat lower than that deduced from 
\J Ej p = 8 X 10^ cm/s, and a Debye lattice specific heat, 
it results Tac < 10~^^ s; Top can be assumed to be of the 
same order of magnitude or even smaller. Since our mea- 
suring frequency is only u) ~ 10^ — 10^, the wt^s factor 
makes the imaginary part negligibly small, and therefore 
the absorption peaks around 100 K cannot be attributed 
to an Akhiezer type mechanism involving the anharmonic 
optic phonon. Still, the real part of 8E contributes to the 
modulus softening with a substantial 5E/E ~ — A(r); 
the continuous curve in Fig. 1 is calculated with the 
above parameters, assuming a temperature independent 
7 = 7. 

Although the mechanisms producing the acoustic ab- 
sorption below 150 K remain unclear, it is worth checking 
whether they can also be the origin of the anomalous soft- 
ening observed below 400 K. The imaginary part of the 
dynamic modulus should be describable in terms of a su- 
perposition of relaxation processes like Eq. (|^), where 
the phenomenological parameter a < 1 describes the 
broadening of the spectrum of relaxation times around 
an average value t, and the relaxation strength is gen- 
erally A oc l/T or nearly constant. In order to repro- 
duce the experimental curves, we need a A (T) strongly 
increasing with temperature, and therefore, we adopt 
A(T) oc 1/ [Tcosh^ (£;/2T)], appropriated for relax- 
ation occurring between states which differ in energy by 
E. The solid lines roughly interpolating the absorption 
data in Fig. 2 are obtained with the following parame- 
ters, T = 5 X 10-" exp (800/T) s, a = 0.2, E/k^ = 230 K 
for the peak at 70-80 K, and 5 x lO^i^ exp (2900/T) s, 
a = 0.4, E/kB = 130 K for the peak at 150 K. We 
do not attribute a particular meaning to these values, 
since it is possible that the two peaks are superimposed 
to a broader maximum, which has not been considered. 
The fit is therefore only indicative, but it shows that the 
anomalous softening below 400 K may be accounted for 
by the same relaxation processes which cause absorption 
at low temperature; in fact, the continuous line in the 
upper panel of Fig. 2 is the real part corresponding to 
the curves in the lower panel and even exceeds the ex- 
perimental softening. The fact that the magnitude of the 
negative step of E' is much larger than the amplitude of 
the E" peaks is a consequence of the non additivity of 
the amplitudes of the elementary peaks, whose maxima 



occur at different temperatures. 

The presence of these intense relaxation processes also 
masks the negative jump and positive change of the 
slope of the elastic moduli due to the superconducting 
transitionp in fact, the latter are of the order of lO^'' 
or less, while the low temperature relaxation processes 
cause a modulus defect of the order of 10"^ at 40 K, and 
they are affected by the transition to the superconducting 
state, as observed in both the real and imaginary parts of 
E. In order to analyze the intrinsic changes of the mod- 
uli at Tc one should be able to subtract the relaxational 
contribution. For this reason, although we can see the 
expected jump and change of slope of E' (T) below Tc, 
we omit the discussion of this issue. 

In conclusion, we measured the complex dynamic elas- 
tic modulus of MgB2 with the principal aim of searching 
for possible lattice instabilities, which might contribute 
to enhancing the electron-phonon coupling without being 
yet revealed by other spectroscopics. The Young's mod- 
ulus indeed presents an anomalous plateau on cooling 
between 400 and 150 K, and hysteresis between cooling 
and heating. The softening, however, can be accounted 
for by intense relaxation processes with maxima between 
50 and 150 K, whose nature is not yet clarified. The 
hysteresis is explainable in terms of generation and anni- 
hilation of geometric defects due to the stresses generated 
by the anisotropic thermal expansion between the grains 
during thermal cycling. The influence of the highly an- 
harmonic vibration modes of the B atoms on the elastic 
and anelastic properties of MgB2 is discussed in detail; 
such anharmonic phonons produce negligible acoustic ab- 
sorption at the low frequencies of the present study, but 
they may substantially contribute to the total softening 
of the elastic modulus with increasing temperature. 
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FIG. 1. Real and imaginary parts of the anelastic spectrum 
of MgB2, measured exciting flexural vibrations at 5 kHz. The 
arrows and numbers indicate the order of the cooling and 
heating runs. The dotted lines in the lower and upper panel 
are a fit to the high temperature relaxation peaks and the cor- 
responding modulus defect, respectively. The dashed line is 
the estimated softening from the anharmonic E2g optic mode 
of the in-plane B vibrations, with the value of the Griineisen 
parameter measured at room temperature. The gray curve 
is the experimental Young's modulus after subtraction of the 
above contributions. 



FIG. 2. Real and imaginary parts of the anelastic spectrum 
of MgB2, measured exciting flexural vibrations at 5.7 kHz and 
73 kHz during the same cooling run. The continuous lines in 
the lower panel are a fit assuming broad thermally activated 
relaxation processes with the parameters given in the text; 
the line in the upper panel is the corresponding contribution 
to the real part of the elastic modulus. The arrows indicate 
the superconducting , at which the curves change the slope. 
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